Monte Carlo Dynamics of driven Flux Lines in Disordered Media 
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We show that the common local Monte Carlo rules used to simulate the motion of driven flux 
lines in disordered media cannot capture the interplay between elasticity and disorder which lies at 
the heart of these systems. We therefore discuss a class of generalized Monte Carlo algorithms where 
an arbitrary number of line elements may move at the same time. We prove that all these dynamical 
rules have the same value of the critical force and possess phase spaces made up of a single ergodic 
component. A variant Monte Carlo algorithm allows to compute the critical force of a sample in a 
single pass through the system. We establish dynamical scaling properties and obtain precise values 
for the critical force, which is finite even for an unbounded distribution of the disorder. Extensions 
to higher dimensions are outlined. 
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In the last few years, the study of elastic manifolds in 
random media has retained much attention. 

These systems appear in a wide range of physical sys- 
tems, ranging from vortices in type-II superconductors 
E|, to charge density waves 0, interfaces in disordered 
magnets ||, and to the problem of directed polymer 
growth gj. The response of elastic manifolds to an exter- 
nal driving force / is highly non-trivial: at temperature 
T = 0, the manifold is completely "pinned" at small 
forces, while it moves with non-zero velocity at forces 
larger than a certain critical force f c . At finite, but small, 
T, a socalled "creep motion" takes place for / <C f c , while 
the motion at / ^> f c is described by viscous flow. Many 
details of this dynamical problem, both at T = and at 
finite temperatures, have yet to be understood fully 

This paper is concerned with an analysis of the dynam- 
ical Monte Carlo method |tJ as applied to lattice models 
of driven elastic manifolds in random media. We argue 
that the common local Monte Carlo rules p|-pl| are in- 
compatible with the Langevin dynamics [Jl2|-|14| which 
defines time evolution in continuum models. We instead 
propose generalized Monte Carlo algorithms where an 
arbitrary number of elements may move at the same 
time. For this class of algorithms, we can establish the 
uniqueness of the critical force and simple connectedness 
of phase space. Furthermore, we devise a method which 
simplifies enormously the calculation of the critical force. 

Our model is sketched in figure [|. We consider a flux 
line x l = {x*}i = o,... ,l moving at times t = 0, 1, 2, . . . on 
a spatial square lattice with a random potential V(i,x) 
with x = 0, 1, . . .. Mainly for convenience (c/ below), we 
also introduce a metric constraint 

\x\ +x -x\\<l (1) 

as well as periodic boundary conditions (x Q = x L ) on the 



flux line. The random potential satisfies 

V(i, x t + M) = V(i + L, Xi ) = V(i, Xi ). (2) 

This condition defines an effective sample of size (L, M). 
While eq.(j2|) is essential for the following, the periodic 
boundary condition for the flux line and the specific 
choice of lattice are inessential details. 

The energy of a flux line x l in presence of an external 
driving force / is given by 

L 

Fix') = J2 V(i, x\) - fx\ + C -{x\ +1 - x<) 2 , (3) 
where c is an elastic constant. 
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FIG. 1. Flux line x = {:E;}i=o,... ,l on a spatial square lat- 
tice with disorder potential V(i, x). Periodic boundary condi- 
tions in i and x are assumed for the lattice and the disorder, 
respectively. The direction of the driving force / is indicated, 
as well as a proposed local Monte Carlo move. The local 
dynamics lead to a trivial critical force in the limit of large 
systems. 

In figure [l], a local Monte Carlo move is indicated. In 
the local Monte Carlo algorithm, the proposed configu- 
ration x differs from the present configuration x only on 
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a random position i. One chooses x~i = X4 ± 1 with equal 
probability. At zero temperature, the move is accepted 
(x t+1 = x) if the energy eq.(||) decreases and if the met- 
ric constraint eq.(|l]) is satisfied. Otherwise, it is rejected 
(x t+1 =x*). 

This rule has been used in past simulations, in spite of 
its very serious shortcomings. Consider, for example, the 
site in figure [l] marked with a circle, (io,Xo). The flux 
line x l shown in figure [l] can only move away from (io, xq) 
if / > V(io,Xo + 1) — V(io,Xo) + c. Even an infinitely 
long flux line (L — ► 00) is thus stopped by a single deep 
pin V(io, Xq) and the motion does not differ qualitatively 
from the one of a point in a disordered potential Jl5],|l6| . 
For an unbounded distribution of V, the critical force 
is always infinite, exposing clearly the pathology of the 
local Monte Carlo algorithm. 

Some authors have therefore used a bounded distri- 
bution, V < V max . However, it is easy to see that in 
this case the local Monte Carlo algorithm does not cor- 
rectly incorporate the disorder, and the flux line's motion 
is trivial: in the limit of large system size LM 3> 1 , the 
flux line will be free to move if / > 2V max , in a way which 
will be similar to the motion of the non-disordered sys- 
tem. This has already been pointed out for the analogous 
case of the random Ising model M . 

We conclude that the description of a driven flux line 
by means of a local Monte Carlo algorithm or its vari- 
ants p"l[ eliminates the very feature which makes the 
problem interesting in the first place, namely the compe- 
tition between elasticity and disorder. This competition 
is preserved in the continuum Langevin dynamics p3|,p|. 

Within the Monte Carlo method, we are thus naturally 
lead to consider generalizations of the model. It can be 
seen easily that the problem just discussed persists even 
if we abandon the metric constraint eq. (Q) . The only re- 
maining option is therefore to abandon the local moves 
in favor of rules which allow global moves. The study of 
global moves in dynamical Monte Carlo is the subject of 
this paper. 



Let us define "Model a" dynamics by a proposed move 
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FIG. 2. Non-local Monte Carlo moves which are considered 
in this paper. Model a: all positions x\ (i = 1, L) may change 
at the same time by a value ±1. Model b: as in a, but the 
motion is either in forward or in backward direction. Model 
c: as in b, but the motion is restricted to single 'fronts'. 



<S* with <r = {61} t 



L such that 



hi p 
8i= { with prob. l-2p ) Vi (Model a) (4) 
-1 p 

At zero temperature, the proposed move is accepted, 
x t+1 — x l + 5, if the resulting configuration both satis- 
fies the metric constraint eq.(jl]) and decreases the string 
energy eq.(^). Note that under Model a dynamics a 
move is proposed with the same probability as its inverse. 
This serves to enforce detailed balance, which allows to 
naturally extend the rule to finite temperatures via the 
Metropolis algorithm. The same can usually not be done 
for cellular automata methods |l3|,[l4|,[l7| . 

A possible second rule ("Model b") chooses at each 
time t with equal probability either to move forward 
(/j* = 1) or backwards (/i f = — 1). The following move is 
then proposed: 



M P 
Si = ( with prob. 

1-p 



Vi (Model b) (5) 



The simulation of dynamical models with such global 
moves may appear hopeless because of the difficulty to 
detect the few energetically favorable choices among the 
exponential number of possibilities in eq.(^) or in eq.(|^). 

To show that the situation is much less desperate, let 
us first define a 'forward front' of length k as a contigu- 
ous set of points i,i+l,....,i + k— 1 which may advance 
together without violating the metric constraints eq.(Q) 
(Si = 5 l+ i ... = S l+k -i = 1 with 5i-i ^ 1, 5 i+k ^ 1). A 
'backward front' is defined analogously. We call 'unsta- 
ble' a front which lowers the energy eq.([|). The moves 
proposed in figure and each consist of two fronts. 
At least one of these must be unstable if the move is to be 
accepted (this is immediately apparent for Model b and 
follows for Model a from an elementary consideration). 
To determine whether a configuration x is unstable, we 
only need to consider the at most 2L(L — 1) + 2 fronts of 
x l rather than the exponential number of moves in eq. (^) 
or eq.(pj). 

Besides Model a and Model b dynamics, it is also pos- 
sible to set up single-front dynamical rules which respect 
detailed balance. These rules (as sketched in figure ^ 
c) can be simulated with less effort than Model a or 
Model b. Even in the latter cases, though, we have devel- 
oped methods which realize eqs (|]) and (^|) while never 
attempting a move forbidden by eq. ([!]). The question 
of which dynamical rule is most satisfactory both from 
physical and computational viewpoints will be dealt with 
elsewhere Q. 

Our main point in the present paper is that a great 
deal of information is available without actually simu- 
lating the dynamic rules. We will show that f c is the 
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same for all models and that the critical flux lines can be 
obtained easily. 

We define, for an arbitrary flux line x a , the 'depinnfng 
force' fd(x a ) as the smallest non- negative / in eq.(||) 
which destabilizes one forward front. Furthermore, we 
define the critical force of a given sample (of size L x M) 
as 

/ c = max/ d (x Q ) (6) 

{x<*} 

where {x a } is the set of all possible flux lines. Notice 
that the definiton of f c or fd(x a ) is model-independent. 
We show in the following that f c in eq.(|^) is an appropri- 
ate definition for all cases as, for a driving force f < f c , 
the system will be pinned in the long-time limit t — > oo. 

To prove the above, we introduce a Variant Monte 
Carlo (VMC) algorithm which, as a byproduct, will allow 
us to actually compute f c with great ease. 

At each time-step t = 0,1,..., the VMC algorithm 
simply moves a single front of minimal length k among 
the unstable forward and backward fronts. The VMC 
method violates detailed balance and is therefore not a 
valid Monte Carlo algorithm. However, each move pos- 
sible within the VMC algorithm is also allowed with all 
the other models considered. 

We have proven the following theorem: if, under VMC 
dynamics at driving force /, a flux line x a is pinned in 
forward direction, it can at most recede towards a config- 
uration x@ [x? < xf Vi), which is itself pinned in forward 
direction. Eventually, we will reach a flux line x 1 which 
is pinned both in forward and in backward directions. 
This flux line x 1 is pinned for all models; if it is pinned 
at / c , we call it a 'critical flux line' x c . 

The theorem allows us to understand that eq.(^|) is 
indeed an appropriate definition for all models: As we 
defined fd{x a ) only with respect to forward motion, one 
might have imagined that a flux line which cannot ad- 
vance at f c could move backwards and then be avoided 
during the subsequent forward evolution. Our theorem 
tells us that such loopholes do not exist: Under VMC 
dynamics, a flux line which can no longer move forward, 
will move backwards and then stop. 

Conversely, we can show that a flux line which can no 
longer move backwards under the VMC dynamics will 
exclusively move forward and then stop. This observa- 
tion simplifies the numerical computations of the VMC 
algorithm. 

Now, we treat the question of how to determine f c and 
x c . There is no a priori guarantee that a generic dynamic 
rule (such as Model a or Model b) will actually encounter 
an x c , when driven at forces / < f c . Simulations in small 
systems, where f c and x c can be obtained by exact enu- 
meration, show that x may pass the sample many mil- 
lions of times without getting pinned. We initially even 
suspected that x c could be dynamically inaccessible from 
part of phase space. 



In this context, we have proven a second theorem: 
Starting at an initial configuration x to with x l ° < x\, the 
VMC algorithm at driving force f < f c can never pass 
x c . In practice, we simply update the driving force by the 
present depinning force / = fd(x ) each time we get stuck 
at a configuration x . In one pass through the system, we 
will have obtained the critical force. The computation of 
f c and of the critical flux line is thus extremely simple. 

Furthermore, the VMC algorithm gives an explicit con- 
struction - for any of the methods - which dynamically 
connects an arbitrary initial state with a critical flux line. 
This proves that all the models in figure || possess a single 
crgodic component. 

We now present our numerical calculations which show 
that f c is finite for a sample of size (L, M) in the limit 
L, M — ► oo with L/M = const. In all our calculations, 
we have used a Gaussian normal distribution for the ran- 
dom potential. 

For finite sizes (L, M), we define the integrated distri- 
bution function Pl,m(/) as the probability that a sample 
of size (L, M) possesses a critical force f c < f. Because 
of the metric constraint eq.(Q), we know that the lateral 
extension of the flux line, max^ j \xi —Xj\, will be at most 
L/2. This can be used to show that for M 3> L/2 

Pl,2mU) = PImU)- (?) 
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FIG. 3. Main figure: Probability Pl,l(/) to have a crit- 
ical force f c < f vs f for samples of size (L, L) vs L for 
L = 16,32,... ,512. The elastic constant is c/2 = 1, and 
a Gaussian normal distribution is used for the disorder po- 
tential. For all large M, the curves Pl,m(/) and P2l,2m(/) 
intersect at the same force, /l,2l which, in the inset, is plotted 
vs 1/L. The extrapolated value of /l,2l in the limit 1/L — *• 
is the critical force f%° of a macroscopic sample. 

We will be interested in the intersection point, /l,2L, 
between the integrated probability distribution for the 
system of size (L, M) and the one of size (2L, 2M): 

Ph,M{fh.2h) = ^2i,2M(/i,2i)- (8) 

In fact, fh,iL will not depend on M for large M, by virtue 
of eq.ln). In our opinion, this observation implies that 
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the natural scaling for our system in the thermodynamic 
limit is M ~ L 1 with 7 = 1, i.e. that we should compare 
the system of size (L, M) with another one, double in 
size both in L and in M. 

We have checked numerically that corrections to the 
scaling relation eq.(^) are already negligible for L ~ M 
(for L > 4) and that intersection points /l,2L indeed do 
not depend on M > L. In figure || we show data for 
PhAf) for i = 16,32, ... , 512. The inset of the figure 
gives the /l,2L as a function of 1/L for all sizes. It is evi- 
dent that Jl,2L extrapolates to a finite value, the critical 
force of the model in the thermodynamic limit, f^°. We 
find f Lf2 L -> !? = 0.91 ± 0.01 for c/2 = 1. We stress 
again that /£° is independent of the aspect ratio L/M. 

For the critical flux line x c , we also studied the size k 
of the minimal unstable front. Naturally, the probability 
distribution of k, p(k), is unbounded, as a finite support 
of p(k) would lead us back to the inconsistencies of the 
local Monte Carlo algorithm. For large k, we find an ex- 
ponential distribution p{k) ~ exp(— k/kt yp ) where kt vp 
depends on the elasticity parameter c, but remains finite 
as L — > 00. Of course, this initial seed of the motion be- 
yond x c may well trigger motions on much larger scales. 
These problems will be studied separately Q . 

Finally, we discuss possible extensions of the work pre- 
sented here. We already indicated that our metric con- 
straint was introduced mainly for convenience and that 
all the developments remain valid. In the absence of the 
constraint, the lateral extension of the flux line may how- 
ever scale as L 7 with 7 > 1. If so, our scaling assumption 
eq.(Q) would have to be modified. We have also extended 
most of our results to higher-dimensional manifolds and 
embedding spaces. There, the only critical issue seems to 
be the complexity of the VMC algorithm, as the number 
of possible fronts can be much larger than in the linear 
flux line. 

In conclusion, we have put the dynamical Monte Carlo 
algorithm for the motion of elastic manifolds in random 
media on a solid footing. We have shown that only 
global-move schemes can capture the subtle interplay be- 
tween elasticity and disorder, which is totally absent from 
the customary local algorithms. Our theorems allowed 
us to compute features universal to all members of this 
class, namely the critical force, as well as properties of 
the critical flux lines. The variant Monte Carlo algo- 
rithm is crucial in that it allows to compute the critical 



force with full rigor even for samples which are several 
orders of magnitude larger than those accessible to exact 
enumeration methods. 

In the future, we think that the extensions to higher 
dimensions will be a very interesting subject for further 
study. Another challenge will be to understand the ac- 
tual time dependence, both at zero and at finite temper- 
atures. It would also be very important to have rigorous 
mathematical proofs concerning the finiteness of the crit- 
ical force and the scaling, in various dimensions. We are 
convinced that the VMC algorithm is sufficiently simple 
to allow such an approach. 
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